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Abstract 

We report on the calculation of the ground-state atomic kinetic energy, E^, and momentum dis- 
tribution of solid Ne by means of the diffusion Monte Carlo method and Aziz HFD-B pair potential. 
This approach is shown to perform notably for this crystal since we obtain very good agreement 
with respect to experimental thermodynamic data. Additionally, we study the structural proper- 
ties of solid Ne at densities near the equilibrium by estimating the radial pair-distribution function, 
Lindemann's ratio and atomic density profile around the positions of the perfect crystalline lat- 
tice. Our value for Ef- at the equilibrium density is 41.51(6) K, which agrees perfectly with the 
recent prediction made by Timms et al., 41(2) K, based on their deep-inelastic neutron scatter- 
ing experiments carried out over the temperature range 4 — 20 K, and also with previous path 
integral Monte Carlo results obtained with the Lennard-Jones and Aziz HFD-C2 atomic pairwise 
interactions. The one-body density function of solid Ne is calculated accurately and found to fit 
perfectly, within statistical uncertainty, to a Gaussian curve. Furthermore, we analyze the degree 
of anharmonicity of solid Ne by calculating some of its microscopic ground-state properties within 
traditional harmonic approaches. We provide insightful comparison to solid 4 He in terms of the 
Debye model, in order to size the relevance of anharmonic effects in Ne. 

PACS numbers: 61.50.Ah,67.80.-s,67.90.+z 
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I. INTRODUCTION 



Noble gases like He, Ne, Xe and Ar, have been intensively studied during the last decades, 
both experimentally and theoretically.-^ Due to their simple electronic closed-shell struc- 
ture, they appear to be affordable many-body systems where to carry out feasible quantum 
computations and test novel methods of calculation. Even though most of them are regarded 
as classical systems, microscopic quantum approaches are required to understand the be- 
havior of the lighter ones, He and Ne, at low temperatures. As it is well-known, He is the 
most representative of the quantum many-body systems. Unique features like Bose-Einstein 
condensation and superfluidity take place in the liquid at few K and recently several ex- 
perimental groups have detected superfluid signal in the solid phase in the mK range,-£ a 
signal that in the homogeneous crystal has been ruled out by accurate theoretical calcu- 
lations.- 1 ^ Moreover, the atomic momentum distribution, n(k), of 4 He differs significantly 
from those of classical systems leading to a non-Gaussian n(k) curve sharply peaked around 
k = 0.-&^ The reasons for those phenomena to happen in helium are the light mass of the 
atoms, bosonic nature of the system and weakness of the interparticle interactions. On the 
other side, Ne has long attracted the interest of condensed-matter scientists since it is an 
intermediate quantum system which provides valuable physical insight when compared to 
other quantum and classical systems. Indeed, the De Boer quantum parameter— defined as 

V mea z 

where m, e and a are the atomic mass, energy scale of the atomic interactions and typical 
interatomic distance of the system, respectively, amounts to 0.54 in Ne (2.50 in 4 He) while in 
Ar and other heavier noble gases, where classical behavior is expected, it drops significantly 
to zero. Essentially, the quantum character of liquid and solid Ne is evidenced on their atomic 
kinetic energies and momentum distributions, which differ appreciably from the predictions 
made by Classical Statistical Mechanics. Accordingly, anharmonic effects in the crystal may 
develop important at low temperatures due to the large zero-point motion of the atoms.— 
In this work, we study solid Ne at zero temperature by means of the diffusion Monte 
Carlo method (DMC)^^ 4 -' 1 ^' and the Aziz HFD-B pair potential.™ Our approach is micro- 
scopic and exact in the sense that the total and partial ground-state energies of the crystal 
may be calculated within statistical uncertainty only. There are burdens of theoretical 
and experimental papers dealing with the thermodynamics and lattice dynamics of solid 



neon, however, numerical results for the atomic kinetic energy are not so abundant. By 
the beginning of the 60's, Bernades 17 and Nosanow et alM were the first in attempting to 
estimate E^ theoretically. They used uncorrelated single-particle wave functions within the 
variational and Hartree approaches, respectively, and arrived at reasonable values not too 
far from present-day calculations; however, the binding energies that they reported were in 
significant disagreement with experimental data. These results made evident the need of 
improved theoretical schemes where to account for the atomic correlations in Ne. Few years 
after Bernades and Nosanow works, Koehler estimated Ej, = 42.6 K by means of the Self 
Consistent Phonon approach (SCP), improving mildly the agreement with experiments.— 

On the experimental side, however, it has not been until the beginning of the 80's, 
with the development of the deep-inelastic neutron scattering technique (DINS), that direct 
measurement of Ef~ in the condensed phases of matter has become accessible. Peek et 
al. performed the first measurements in solid Ne, covering the temperature interval 4.5 — 
26.5 K.— The authors of the first study reported E^ = 49.1 ± 2.8 K for the ground-state 
kinetic energy and, because of the large discrepancies with respect to calculations based on 
harmonic models, they suggested substantial anharmonic effects in solid Ne. 

Reassuringly, few years after Peek et a/.'s measurements,— theoretical estimations by 
Asger and Usmani,— who used a perturbational approach based on a Wigner-Kirkwood high- 
temperature expansion with the Lennard- Jones (L-J) and Aziz HFD-C2 24 pair potentials, 
amounted to E^ ~ 49 K at temperatures near 10 K. Regardless, previous to Asger et 
a/.'s results,— Cuccoli and co-workers^ 2 - arrived at kinetic energies ~ 7 K below Peek's 
results, based on the full quantum approach path integral Monte Carlo (PIMC) and the 
L-J interaction. The authors of this work suggested that their disagreement with Peek's 
results could be in part due to the oversimplification of the atomic interactions made by the 
adopted potential. Subsequently, Timms et al— performed a series of new low-temperature 
DINS experiments in solid Ne at high momentum tranfers with an improved experimental 
set-up. They found very good agreement with Cuccoli et al— and also with Ceperley and 
Boninsegni, 23 who performed an exhaustive PIMC study of the crystal at low temperatures 
using both L-J and HFD-C2 pair potentials. In addition, a recent theoretical study by 
Neumann and Zoppi, in which computational techniques and interatomic potentials similar 
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are used, comes to reinforce the accuracy of Timms' data.— Very recently, 



to those of Ref. 

Timms et al?- have reported new additional DINS measurements in solid Ne performed 
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within the temperature range 4 — 20 K. By doing this, they complement their previous 
results and provide a truster way to infer the value of Ek in the ground state, which by 
means of extrapolation of the excess kinetic energy turns out to be 41(2) K. 

In the present work, we report quantum Monte Carlo results of the equation of state 
and other thermodynamic properties of solid Ne over a range of densities near equilibrium 
(—1.2 < P < 6 Kbar), and find overall excellent agreement with experimental data. Struc- 
tural properties of the crystal, namely the radial pair-distribution function, g(r), atomic 
density profile around the positions of the perfect crystalline lattice (sites) and Lindemann 
ratio, are also provided. Remarkably, we estimate accurately the atomic kinetic energy 
of the crystal at its equilibrium density by means of the pure estimator technique within 
DMC— i2&2£ Our result, Ek = 41.51(6) K, is in very good agreement with the recent predic- 
tion of Timms et al. 2 ® We have also calculated the ground-state atomic momentum distri- 
bution n(k) of solid Ne and it is found to fit perfectly to a Gaussian within the statistical 
uncertainty. 

Additionally, we have analyzed the degree of anharmonicity of solid Ne in its ground 
state. With this aim, we have computed the atomic kinetic energy and mean squared 
displacement within the Self Consistent Average Phonon (SCAP) approach,— 1 ^ which is 
a simplified version of the Self Consistent Phonon method^ that has proved successful in 
reproducing a deal of thermodynamic properties of rare gas solids. We find the SCAP results 
are not in full agreement with the quantum DMC ones, thus revealing this approach might 
not allow for an accurate description of Ne at the microscopic level. In a further step, we 
devise an harmonic model based on the HFD-B potential in which the interaction between 
particles depend on their relative distances, equilibrium positions and the force constant field 
(second derivatives of the potential energy evaluated in the perfect crystal configuration). 
By using DMC, we calculate the total and kinetic energies associated to this model and 
find significant discrepancies with respect to the full HFD-B results. According to these 
outcomes, solid Ne may be regarded as a moderate anharmonic crystal since, contrarily 
to what is observed in solid 4 He, its n(k) does not deviate appreciably from the Gaussian 
pattern. 

The remainder of this article is as follows. In Sec. [Til we describe the computational 
techniques and models that have been used on this study. Next, in Sec. IHH we present our 
results and compare to previous experimental and theoretical data. In Sec. IIV1 we finalize 
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by summarizing the main conclusions and giving some general remarks. 



II. TECHNIQUES AND MODEL 



A. Diffusion Monte Carlo 



DMC is a zero-temperature method which provides the exact ground-state energy of the 
many-boson interacting systems within some statistical errors.— This technique is based 
on a short-time approximation for the Green's function corresponding to the imaginary 
time-dependent Schrodinger equation, which is solved up to a certain order of accuracy 
within an infinitesimal interval At. Despite this method is algorithmically simpler than 
domain Green's function Monte Carlo,— >^ it presents some (Ar) n bias coming from the 
factorization of the imaginary time propagator e~T H . Nevertheless, our implementation of 
DMC is quadratic,— hence the control of the time-step bias is efficiently controlled since the 
required At — > extrapolation is nearly eliminated by choosing a sufficiently small time 
step. The Hamiltonian H, describing our system is 

- 2 n N 
H--^£V? + 5>M, (2, 

1Ne 1=1 Kj 

where m-Ne is the mass of a Ne atom, the distance between atoms composing an i,j pair 
and V^Tij) the interatomic interaction that we have chosen as the Aziz HFD-B potential.— 
The corresponding Schrodinger equation in imaginary time (it = r), 

- = (H - E) *(R,r) (3) 

with E an arbitrary constant, can be formally solved by expanding the solution \1/(R, r) 
in the basis set of the energy eigenfunctions {$ n }- It turns out that \1/(R, r) tends to the 
ground-state wave function $ of the system for an infinite imaginary time as well as the 
expected value of the Hamiltonian tends to the ground-state value E . The hermiticity of 
the Hamiltonian guarantees the equality 

($o|H|$ ) Ngjr) /m m 
($ |$o) ($o|V»t) 

where ipT is a convenient trial wave function which depends on the atomic coordinates of 
the system R = {ri,r 2 , ...jTn}. Consequently, the ground-state energy of the system can 



be computed by calculating the integral 

(E) DMC = lim / E L (R) / (R, r) dR , (5) 

where / (R, r) = \& (R, r) (R), and £x (R) is the local energy defined as Ex(R) = 
H/0y (R) /^t (R)- The introduction of (R) in / (R, r) is known as importance sampling 
and it certainly improves the way in which integral (111) is computed (for instance, by imposing 
ipT (R) — when is smaller than the core distance of the interatomic interaction). 

In this work, all the operators diagonal in real-space which do not commute with the 
Hamiltonian, that is [H, O] 7^ 0, have been sampled with the pure estimator technique.— i22i22 
With this method, essentially, the possible bias induced by ipT in the mixed estimator 
(^&o\0\tpTj are removed by proper weighting of the configurations along the simulation. 



B. Trial wave function and pair potential 

We have modeled solid Ne by assuming point-like atoms interacting via a radial pair-wise 
potential and with equilibrium positions distributed according to the fee structure. Neon is 
observed to remain stable in the fee structure up to pressures of 1100 Kbar and at ambient 
temperature,— therefore, no other configuration apart from this has been considered in the 
present study. The potential chosen for the interatomic interactions is the semi-empirical 
Aziz HFD-B one,— which has proved excellent in reproducing some of the macroscopic and 
microscopic properties of Ne over a wide range of temperature and pressure, 36 and appears 
to be more realistic than the Aziz HFD-C2^ and Lennard- Jones (L-J) models at short 
distances. Explicitly, 

V(r) = eB(x) , (6) 

where 

Q(x) = Ae W (-ax + ^)~F(x)^ + ^ + ^) , 

and 

F(x) = ( <*» [- - * < D 
[l x > D . 

The value of the parameters of the potential are A = 895717.95 , a = 13.86434671, D = 1.36, 

r m = 3.091 A , p = -0.12993822 , e = 42.25 K , C 6 = 1.21317545 , C 8 = 0.53222749 and 

C10 = 0.24570703 , with x = r/r m . It is known that, upon high pressure the introduction 
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of additional terms in the effective atomic potentials of rare gases are required to account 
for many-body effects taking place on them; for instance, in solid Ar this limit is posed 
around 50 Kbar.— ^ This circumstance, however, does not affect the reliability of the results 
that we are to present in short, since the pressure range involved in our simulations is 
-1.2 < P < 6 Kbar. 

Regarding the trial wave function chosen for importance sampling, ipT, we have adopted 
the extensively used and tested Nosanow-Jastrow model,— ^241 

TV N 

i, T (n, r 2 , r N ) = llHnj) ngidri - Ri|) , (7) 

i^j i=i 

with f2(V) = e - ^) a nd gi(r) = e~^ ar2 . The best parameter values are a = 6.5 A -2 , 
b = 4.0 A and c = 5.0 , optimized using the variational Monte Carlo method. Their 
dependence with the pressure is small and therefore neglected for its use on the DMC 
simulations. The first factor in ip? accounts for the correlations between particles induced 
by the interactions, while the second enforces the atomic ordering within the system by 
attaching each particle to one site of the perfect lattice through a Gaussian function. The 
indistinguishability of the Ne atoms has been neglected throughout this work since the 
Nosanow-Jastrow model is not symmetric under the exchange of particles. This choice is 
fairly justified since quantum effects derived from a correct symmetrization are not expected 
to play any significant role in the solid properties calculated in this work. In fact, the same 
conclusion for the same quantities also holds for solid 4 He, a solid with a larger quantum 
behavior. The parameters of the simulation, namely the number of particles per box, time 
step and target walker population (that is, the mean number of walkers along the simulation), 
have been chosen in order to ensure the correct asymptotic behavior; their respective values 
are: N = 256 , Ar = 2.7 • 10~ 4 K" 1 and n w = 260 . 

At each density, finite size effects have been corrected by including the tails of the kinetic 
and potential energies into the total energy, both estimated assuming g(r) = 1 beyond half 
the length of the simulation box. This assumption could be too crude for solids (see Fig. [2]) 
and therefore we have checked the reliability of this approximation in our system. To this 
end, we have carried out some simulations with 500 atoms and compared the energetic and 
structural results with the ones obtained for a box of 256 particles at the same density. For 
instance, at a density p = 0.045 A" 3 the energies are E/N = -238.88(4) K and -238.69(8) K 
and the Lindemann ratios 7N e = 0.077(1) and 0.079(3) for N = 256 and 500 particles, 
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respectively. The differences observed are therefore not significant within our statistical 
uncertainty and the size corrections are reasonably included. 



III. RESULTS 

A. Thermodynamic properties 

In Fig. [1] (left), we show our results for the total atomic energy of solid Ne at zero 
temperature. The solid line on it corresponds to the polynomial curve, e(p) = E (p) /N, 

which has been fitted to the DMC energies reported in Table [J (solid points in the fig- 
ure). The values of the parameters of the best fit are a = 938(3) K, b = 871(20) K, 
e = —239.21(3) K and p = 0.04582(2) A -3 , where e and p are the equilibrium energy per 
particle and density, respectively. The agreement between our results and experiments is rea- 
sonably good for the energy e e xpt = -232(1) K and the density p e xpt = 0.044976(3) A" 3 .-^ 
Once e(p) is known, it is straightforward to deduce the pressure, P{p) (see Fig. [TJ right), 
and compressibility, ft(p), of the system at any density through the relations 



P(P) = P 



2 de(p) 



dp 



= (9) 

The compressibility at the equilibrium density obtained so is = 0.084(4) Kbar -1 , which 
compares excellently to the experimental value K,Q Xpt = 0.089(2) Kbar -1 .— 

An interesting magnitude in the study of condensed phase systems is the spinodal density, 
ps, which is the thermodynamical limit for the system to remain in homogeneous phase. 
At this density, the relation dP/dp = is fulfilled, which is equivalent to require infinite 
compressibility or zero speed of sound in the system. Our prediction for pg is 0.03575(5) A~ 3 , 



which corresponds to a pressure P(ps) = —1.102(4) Kbar. In Ref. [44|, Herrero presents a 
comprehensive study of solid Ne at negative pressures by means of the PIMC method. 
The author modelizes the interatomic interactions with the L-J potential and estimates 
the pressure at the spinodal density and zero temperature by means of a linear fit to the 
squared bulk modulus with respect to pressure; he obtains P(ps) PIMC = —0.91 Kbar and 
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FIG. 1: Left: Energy versus density for solid Ne at zero temperature. The solid line corresponds 
to the polynomial curve of Eq. [8] and the dots to the calculated DMC total energies per particle, 
errors bars are smaller than the size of the symbols. Right: Equation of state of solid Ne at zero 
temperature computed with DMC and the Aziz HFD-B potential. 



(A" 3 ) 


E/N 


E p /N 


E k /N 


0.040 


-225.84(4) 


-256.26(8) 


30.34(8 


0.043 


-235.83(4) 


-272.04(8) 


36.23(8 


0.044 


-237.88(4) 


-276.50(8) 


38.57(8 


0.045 


-238.88(4) 


-279.43(8) 


40.61(8 


0.047 


-238.55(4) 


-283.76(8) 


45.17(8 


0.050 


-230.76(4) 


-282.81(8) 


52.17(8 


0.053 


-212.83(4) 


-272.15(8) 


59.31(8 


0.056 


-183.20(4) 


-249.61(8) 


66.41(8 



TABLE I: Total, potential and kinetic energies per particle of solid Ne at absolute zero as computed 
with DMC and the pure estimator technique. Energies are in units of K. 

Pg IMC = 0.0356 A~ 3 . The disagreement between this and our value for P(ps) can be 
explained in terms of the adopted interatomic potential, since small differences in the total 
energies may develop large within successive derivatives. 
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B. Structural properties 



We have explored several structural properties of solid Ne. In Fig. [2j we plot the averaged 
radial pair-distribution function, g(r), which is proportional to the probability of finding a 
particle at a certain distance r from another. According to what is expected in crystals, 
g(r) emerges peaked with maxima corresponding to the distances between successive shells 
of atoms within the perfect lattice, though the peaks broaden with respect to the profiles 
which are obtained in classical solids. 

A characteristic parameter in the study of quantum solids is the Lindemann's ratio, 7, 
which is defined as the ratio between the squared root of the mean squared displacement, 
(u 2 ), and the distance between first nearest neighbours in the perfect crystalline lattice. 
Our estimation of the Lindemann's ratio at the equilibrium density (pure estimation) is 
7Ne = 0.088(2) , which is significantly smaller than in 4 He (~ 0.26) and H 2 (~ 0.18), but 
still larger than in classical solids at finite temperature and far from melting (~ 0.03). The 
corresponding mean squared displacement, (u 2 e ), amounts to 0.077(1) A 2 . In Table HU we 
quote the value of 7N e at several densities out of the equilibrium. As it is observed therein, 
the general trend of 7N e is to reduce when the density is increased; this behavior is easily 
understood in terms of gain of cohesion energy, which must balance with the increasing of 
kinetic energy of the system arising from atomic localization. 

Aimed to characterize the spatial distribution of the atoms around the equilibrium po- 
sitions in solid Ne, we have calculated the atomic density profile function (averaged for all 
directions), /i(r), and kurtosis, (q . The averaged atomic density profile function, /i(r), yields 
the probability of finding a particle at a distance within the interval (r, r + dr) from any 
arbitrary site of the lattice. According to this definition, the mean squared displacement, 
(u 2 ), can be obtained as 



In Fig. [3l we plot /i(r) at the equilibrium density (dots), together with the Gaussian curve 
that we have adjusted to it (solid line). To check the reliability of this fit, we have assumed 
the Gaussian curve in Eq. ( flOi) . instead of /i(r), and then recalculated (u 2 ). Proceeding so, 
we obtain 0.079(1) A 2 which agrees perfectly with the direct calculation 0.077(1) A 2 . Next, 
we compute (q in several directions of the cubic cell so as to discern whether the atoms 




(10) 
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P (A 3 ) 7Ne C(100) C(OIO) 



0.040 0.099(2) 0.017(12) 0.012(14) 

0.043 0.092(2) 0.000(8) -0.001(8) 

0.044 0.091(2) 0.000(10) 0.000(10) 

0.045 0.087(2) -0.006(7) -0.014(7) 

0.047 0.086(2) 0.000(20) 0.000(10) 

0.050 0.083(2) 0.000(10) -0.010(10) 



TABLE II: Lindemman's ratio, 7Ne 5 an d kurtosis, £q of solid Ne at different densities close to 
equilibrium. 

distribute isotropically in average or not around the sites. The kurtosis is defined as 

\ u (ijk)) 

where viujk) are the projections of the position vectors which relate each lattice site to its 
nearest particle along the {ijk) direction (Cartesian basis). As it is well-known, if the atomic 
density distribution over the equilibrium positions is of Gaussian type the kurtosis is null. 
In the case of solid Ne, we have obtained C(ioo) — 0.0078(63) and C(oio) — 0.0062(59), which 
indeed might be regarded as values compatible to zero. Additional results for (q obtained 
with the pure estimator technique are quoted in Table [III 



C. Kinetic energy and momentum distribution 

In Table HI we summarize the value of the atomic ground-state kinetic and potential 
energies of solid Ne near equilibrium (P ~ 0). All the E p and Ek results have been computed 
within the pure estimator technique and DMC, thus any possible errors associated to them 
are of statistical kind or stem from the modelization of the interatomic interactions. In 
particular, we have estimated E/, = 41.51(6) K at the equilibrium density. In Fig. HJ we 
plot the values of the excess kinetic energy of solid Ne, defined as E exc — E^ — (3/2) T, as 
measured by Peek— and Timms et al.— within the temperature range 4 — 20 K. Therein, we 
also include estimations of E exc as obtained with PIMC over the same T-interval, together 
with our ground-state result which is located at the ordinate axis. By performing linear 
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FIG. 2: Averaged radial pair-distribution function, g(r), of solid Ne at zero temperature and the 
equilibrium density. 

fits to the excess kinetic energy, it is shown that our ground-state prediction is in very 
good agreement with Timms's measurements^ and the PIMC estimations,— whereas not 
so with Peek's results.— The causes for this disagreement may be explained, as it has been 
suggested elsewhere,— in terms of systematic experimental errors, since the temperature 
dependence of E exc obtained by Peek and co-workers appears to coincide with Timms' results. 
A likely explanation can rely on the range of neutron momentum transfers involved in 
those first DINS experiments, about two orders of magnitude less intense than in posterior 
measurements, which might not be sufficiently large so as to reach the high Q-regime required 
for the impulse approximation of the dynamic structure factor to be valid.— 

Another physically rich quantity in the study of quantum liquids and solids is the one- 
body density matrix, g(r,r'), which is defined as 



where ^(r') and -^(r) are, respectively, the field operators which destroy a particle from 
position r' and create one at position r and $o is the ground-state wave function. In 
boson systems the asymptote lim^oo g(r) provides the condensate fraction of the associated 




(12) 
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FIG. 3: Atomic averaged density profile, n(r), of solid Ne at zero temperature and the equilibrium 
density. 

homogeneous system no- The Fourier transform of g(r) is directly the atomic momentum 
distribution, 

n(k) = p J dv e ikr g{r) . (13) 

In the Quantum Monte Carlo formalism, the one-body density function can be estimated by 
averaging the coordinate operator A(r, r 1; r N ) = i/jt( t i + r ; r 2, •••> r N)/V'T( r i, r 2, r N) 
within customary DMC (known as mixed estimation, g m i X {r) = (A(r)) dmc)-~ However, a 
more accurate evaluation of g(r), known as extrapolated estimation, is given by the expres- 
sion 

g(r) = 2g mix (r) - g V Mc(r) , (14) 

where gvMc{ r ) results from averaging A(r, rj., tn) within variational Monte Carlo. In 
Fig. [5], we plot our results for g(r) as given by Eq. (TT^j) . In the same figure, we also enclose 
the Gaussian curve, G(r) = e~ br2 (given that £>(0) = 1), which best fits to our calculations, 
with an optimal parameter value b = 5.743(36) A -2 . In order to test the quality of this fit 
(which in the reduced chi-squared test gives the value 0.99), we have calculated the atomic 
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FIG. 4: Excess atomic kinetic energy of solid Ne at low temperatures. Experimental data of Ref. \2G 
are represented by A , measurements of Ref. 



20 by • , PIMC estimations of Ref. 



23j by ▲ and our 



ground state estimation by T (in the ordinate axis). The lines in the plot correspond to linear fits 



to the experimental data of Refs. 



20 



26 . 



kinetic energy of solid Ne through the formula 

h 2 



2m 



V 2 g(r) 



No 



(15) 



r=0 



but assuming G(r) instead of g(r). In fact, it may be shown that Eq. (1151) derives from the 
kinetic-energy sum rule 



Eh 



h 2 



dk k 2 n(k) 



(16) 



2m Nc (2tt) 3 p 

Proceeding so, we have obtained E^ = 41.43(26) K, which fully agrees with the direct 
estimation 41.51(6) K. This finding allows us to conclude that g(r) in solid Ne at T = can 
be well considered Gaussian-shaped at all effects. 

We have also computed the atomic momentum distribution of solid Ne by taking the 
Fourier transform of g(r) over a set of k- vector points, as expressed in Eq. ( jl~3l) . In Fig. El 
we plot the results of these calculations (dots) and additionally the Fourier transform of the 
aforegiven Gaussian fit to g(r) (solid line with width signalizing the associated uncertainty). 
Obvioulsy, once g{r) has proved Gaussian, n{k) turns out to be of the same kind. 
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0.2 0.4 0.6 0.8 1 1.2 

r(A) 

FIG. 5: One-body density matrix of solid Ne at the equilibrium density. The solid line in the figure 
corresponds to the Gaussian curve that we have fitted to the results. 

D. Degree of anharmonicity 

The Self Consistent Phonon approach^ (SCP) has proved very accurate in characteriz- 
ing solids in the middle way between classical and quantum behavior. Very essentially, this 
theory makes the assumption of particles coupled harmonically with frequencies and modes 
depending on the crystal symmetry and lattice parameter and which are determined through 
a self-consistent procedure. A simplified version of this method is the Self Consistent Aver- 
age Phonon approach (SCAP),— £L which adopts the expressions of SCP but replacing the 
summation over the different vibrational frequencies by an averaged one, namely the Ein- 
stein frequency, Qq. Despite this crude simplification, the agreement between measurements 
and SCAP results for rare gase solids like Ar, Kr and also Ne is excellent in what con- 
cerns thermodynamic properties (isothermal bulk modulus, specific heat, etc.).— However, 
we want to know to what extent harmonic assumptions in solid Ne are accurate enough for 
deriving microscopic properties of its ground state. With this aim, we have calculated the 
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FIG. 6: (Color online) Momentum distribution of solid Ne at the equilibrium density (green dots 
and bars). The solid line in the figure corresponds to the Fourier transform of the Gaussian curve 
previously fitted to g(r) (the width of the line represents the uncertainty of the fit). 

atomic kinetic energy and mean squared displacement within SCAP through the formulas 
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n 2 

"0 



2 r 



(17) 
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and 
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3h 



2m Ne Q 

1 m Ne (u 2 )( SCAp )^ = l m 



o • 



(18) 



2 , u 4 

where V^(r) is the Aziz HFD-B pair potential. We first compute the exact value of Qo with 
the pure estimator technique within the DMC approach, and then calculate the value of 
expressions (ITHjl . The results that we have obtained are, HQ = 62.04(1) K, (ulf c ) SCAP = 
0.058(3) A 2 and (E k ) SCAP = 46.5(1) K, which disagree noticeably from the DMC values 
(u^ e ) = 0.077(1) A 2 and (Ep.) = 41.51(6) K. This outcome reveals that crude simplifications 
made on the vibrational properties of solid Ne may lead to important inaccurracies on the 
quantum description of such crystal. 
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In a further step, we have devised an harmonic model^ in which the interaction between 
particles is pairwise and reads 

Vt 1 ^) = V 2 {r m ) + \{v*- Uj? {^-) K - u,) , (19) 

V ij ij / rij=ro,ij 

where Vi(r) is the Aziz HFD-B interaction, Uj is defined as r— Rj, and the terms V2(ro,ij) and 
( dr~lF~ i m ^ e right side of Eq. ( fT9l) are evaluated, only once, for the atoms in the 

V r y r !j / rij=ro y 

perfect crystal configuration (r .jj = |Rj — Rj|) . This approach is equivalent to assume the 
pair of atoms % and j coupled through an harmonic spring of constant equal to the second 
derivative of V^(r) evaluated at the equilibrium distance r ^j . Within DMC and with 
the pure estimator technique, we have computed the exact ground-state total and kinetic 
energies of this model, arriving at the values, eQ arm = —251.35(4) K and E\ &ria = 35.1(3) K, 
which differ notably from the results obtained with the full Aziz HFD-B interaction. 

The relative failure of the previous approximations allow us to conclude that traditional 
harmonic approximations in solid Ne are not adequate for an accurate evaluation of its 
microscopic properties. Aimed to yield a rough estimation of the degree of anharmonicity of 
solid Ne, and to finalize with this section, we now compare solid Ne with solid 4 He, the most 
anharmonic among all the crystals, by invoking the Debye model. In the Debye approach 
for solids, particles are assumed as non-interacting quantum harmonic oscillators which 
vibrate with frequencies within a spectrum that is top-bounded by the Debye frequency, 
uj d . Consequently, the atomic kinetic energy is expressed as Eg = (9/16) 0£>, where 9^ is 
the Debye temperature and is equal to Hujo- It is readily shown that Qd = 9h 2 /4m(u 2 ) , 
which in the case of solid Ne at equilibrium turns out to be 70.3(9) K (here, we have used 
the value (ug, e ) = 0.077(1) A 2 ), which in turn leads to Eg = 39.5(5) K. Next, we define the 
dimensionless parameter, T = 1.0 — [Eg /EjA, which in fact vanishes for the case of a pure 
harmonic solid (Debye model) and it progressively increases towards unity as anharmonic 
effects develop larger. For solid Ne and 4 He at their respective zero-temperature equilibrium 
volumes, we assess the values T-^ e = 0.05 and Tne = 0.44, where for helium we have used the 



data found in Ref. |47|. By comparing these two figures, one could claim that anharmonic 
effects in solid Ne are about one order of magnitude less substantial than in 4 He. 
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IV. DISCUSSION AND CONCLUSIONS 



In this work, we report the calculation of the ground-state atomic kinetic energy, one- 
body density matrix and momentum distribution of solid Ne by means of the DMC method 
and the realistic Aziz HFD-B pair-potential. Our approach is proved to perform notably 
for this crystal, as it is shown by the very good overall agreement obtained with respect 
to thermodynamic experimental data. Our value for the atomic kinetic energy of solid Ne 



at the equilibrium volume, E\. ■ 
experimental data found in Refs. 



41.51(6) K, is in accordance with the low-temperature 
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26l and also with previous PIMC calculations performed 



with the L-J and Aziz HFD-C2 pairwise interactions.— 1 ^ 1 ^ However, our result does not 
agree with the results obtained by Peeks and co-workers (previous to Timms' work) based 
also on deep-inelastic neutron scattering measurements. We have calculated the one-body 
density function of solid Ne and shown that it perfectly fits to a Gaussian curve. Conse- 
quently, the atomic momentum distribution, which is evaluated by performing the Fourier 
transform of g(r), is of the same kind. Interestingly, Withers and Glyde^ have shown very 
recently by means of simple models that the deviation of n(k) from a Gaussian pattern 
in quantum solids may arise by effect of anharmonicity and/or the introduction of atomic 
exchanges. We have checked that anharmonic effects in the ground-state of solid Ne are rel- 
evant by calculating some of its microscopic properties within traditional harmonic schemes 
and quoting significant discrepancies with respect to the full quantum results. It is noted 
that we have not attempted to include atomic exchange effects in the present work since 
a priori and very reasonably, these are not expected to play any substantial role in the 
ground-state of solid Ne (contrarily to what may occur in 4 He, for instance). Even so, we 
do not appreciate, within the statistical uncertainty, any deviation from a Gaussian pattern 
in the g(r), or equivalently n(k), of solid Ne, therefore, the degree of anharmonicity of Ne 
at zero temperature may be regarded as fairly moderate. 
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